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Abstract 

■ The electrical conductivity in the hot phase of the QCD plasma is extracted from a quenched 

o : 

CN . lattice measurement of the Euclidean time vector correlator for 1.5 < T/Tc < 3. The spectral 

^ I density in the vicinity of the origin is examined using a method specially adapted to this region, 
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and a peak at small energies is seen. The continuum limit of the electrical conductivity, and the 



. closely related soft photon emissivity of the QCD plasma, are then extracted from a fit to the 

> , 

. Fourier transform of the temporal vector correlator. 
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The soft photon production rate from the plasma phase of hadronic matter is of impor- 
tance to searches for the QCD phase transition Consequently, there has been a long 
history of attempts at perturbative computations of this rate [2|. The first lattice prediction 
of dilepton (off-shell photon) rates was performed a while back ||3||1. Recently the leading 
order computation of the photon production rate was completed J]- The Kubo formula 
relates the soft limit of this rate to the DC electrical conductivity of the QCD plasma, a. To 
leading log accuracy in the gauge coupling, g = v^vraj, one has a oc aT / g'^ log , where 



a is the fine structure constant. T 
in the leading-log approximation 



J. Ill 



le proportionality constant has been computed recently 
Here we report the first computation of a and the soft 



photon emissivity from a quenched lattice computation in a region of temperature where g 



is large and the weak-coupling approach fails ^|. Our methods can also be applied to other 
transport problems. 

The photon emissivity at temperature T is related to the imaginary part of the retarded 
photon propagator, i.e., the spectral density, p^*^^, for the electromagnetic current correlator, 
through the relation 

In this work we shall take = p = 0, and hence obtain the soft photon production rate. 
Since the EM Ward identity gives p°°(c<j,0,T) = 0, this soft limit is related to transport 
properties of the QCD plasma through the Kubo formula. 



(2) 



uj=0 

where the sum is over spatial polarisations. A lattice determination of this rate proceeds 
from the spectral representation for Euclidean current correlators — 

GEM{t,p,T) = —K{uj,t,T)pEM{uj,p,T), (3) 
Jo Zn 

where the integral kernel K = exp{u!t)ng{u!,T) + exp{—ujt)[l + rig^UjT)]. Gem is the 
product of the vector correlator summed over all polarisations, Gv, and the EM vertex 
factor Gem = 4:naJ2f^}, where e/ is the charge of a quark of flavour /. On discretising 
the integral it becomes clear that the extraction of Pem from the lattice computation of 
Gem is akin to a linear least squares problem. The complication is that the (potentially 
infinite) number of parameters to be fitted exceeds the number of data points (which is 
half the number of lattice sites in the time direction, A^^^). The solution is to constrain the 



function Pem through an informed guess 
The Maximum Entropy Method (MEM) 



IG 



spectral function has been used in the past 



and use a Bayesian method to extract it 



ll| along with a free-field theory model of the 

'or C(j/T > 4 is fully 
For that reason 



'Sj. The hard dilepton rate 
under control, with lattice and perturbation theory in good agreement 
we concentrate here on the electrical conductivity and the soft photon rate. 

Correlators were investigated at T = 1.5Tc, 2Tc and 3Tc in quenched QCD. The tem- 
perature ran ge i s realistic for heavy-ion collisions. However, g > 1 m this entire range of 



temperature and is therefore ineffective in the separation of length scales upon which 
weak-coupling approaches depend. In order to make continuum extrapolations, the com- 
putations were performed on a sequence of lattice spacings, a = 1/8T, 1/lOT, 1/12T and 
1/14T, (i.e., Nt = 8, 10, 12 and 14). Quark mass effects were controlled by working with 
staggered quarks of masses m/Tc = 0.03 and 0.1. Details of the runs, statistics, and the 
generation of configurations for A^^^ < 14 are described in For these lattice spacings 
the computations were performed on two different spatial volumes in order to control finite 
volume effects. For Nt = 14 we have added runs on 14 x 30^ lattices for T = 1.5Tc and 2Tc, 
and on 14 x 44^ lattices for T = 3Tc, generating 50 configurations separated by 500 sweeps 
each. We have measured vector correlators with two degenerate flavours of quarks. It has 
been demonstrated recently that in this limit the charged and uncharged vector correlators 
are identical jisl . 

Small but statistically significant differences between the lattice results and ideal quark 
gas predictions for Gem are observed at all temperatures, lattice spacings, quark masses and 
volumes investigated. In any lattice computation, we expect the high frequency part of Pem 



to contain lattice artifacts. Moreover, physics at momenta of order 1/a is perturbative [14 1 
and not of interest in the present context. We remove this physics by taking the difference 
between the Euclidean temporal propagators in QCD and an ideal quark gas (free field 
theory) on the same lattice — 

AGEM{uj)=GrM''{u^)-GT.f{u^), (4) 

A subtraction is needed to remove the cj^ divergent pieces from the dispersion relations . 
Other benefits accruing from this are discussed later. 

We have estimated the spectral density by two classes of methods. The first class of 
general techniques consist of discretising the integral in eq. d^) into N^^ energy bins and 
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rewriting the equation in the form Gem = Kp^M where K is now an Nt x A^^ matrix, Gem 
the data vector of length Nt and Pem represents a vector of length A^^^ [16]. For A^^^ > Nt 
the solution is non-unique. Additional constraints, called priors, must then be imposed to 
determine them The extraction of the spectral density is performed in the context of 
Bayesian parameter extraction. Given the data on G, the probability distribution function 
for p can be written using Bayes' formula 

P{p\G) = P{G\p)P{p)/P{G), (5) 

where P{A\B) denotes the conditional probability of A given B. The probability distribution 
function P{p) contains the prior information on p that is needed for the analysis. Writing 
P{p\G) = exp[—F{p)], the maximum likelihood analysis of the probability reduces to the 
problem of minimizing F. Since P{G) is independent of p, this problem can be formulated 
as a minimisation of the function 

F(p) = (G - Kpfj:-\G - Kp) + /5t/(p), (6) 

where the first term is the logarithm of P{G\p) and P{p) = exp[—f3U]. The superscript 
T denotes a transpose, S is the covariance matrix of the data, and (3 is a non-negative 
parameter whose choice is specified later. 

The MEM technique consists of choosing some vector p^^,^ and defining U{pem) = 
PsMl^'^&iPEM/ pVm) ~ 1]' where the sum is over components of the vectors. In previ- 
ous works the prior p^^,^ has been chosen to be the vector spectral function in an ideal quark 
gas p. Another whole class of techniques is obtained by choosing U{pem) = [^{Pem — P^em)]'^ 
where £ is a non-singular matrix. The choices C = 1, D and (where D is a discreti- 
sation of the derivative) have been suggested in the literature. £ = 1 is the model that 
Apem = Pem — P%M = except as forced by the data, C = D makes the a priori choice that 
Apem is constant and C = is the prior choice of smooth Apem 13 • 

Such regulators have the added advantage that minimisation of the function F in eq. ^ 
yields the linear problem — 

K^E~'K + Pem = K^^-'Gem + (3C^Cpl,,. (7) 

Since C is positive definite, it is clear that the term in j3 on the left hand side regulates 
the problem, by adding a term to K^T,~^K which makes the sum invertible. Since, for a 




FIG. 1: Bayesian fits to AGv{t) at T = 2Tc for m/Tc = 0.03 and Nt = 8 (circles) 10 (squares) and 
12 (pentagons). The fits were made with N^^ = 16 and < lo < AttT, choosing £ = 1. Changes in 
the fits due to variations in these algorithmic quantities are indistinguishable on the scale of this 
figure. 

well-determined parameter fitting problem, the value of F is the value, we choose a value 
of f3 at which F = Nt at the minimum of F{f3), i.e., at the maximum a posteriori probability. 

Considering the Bayesian problem as a field theory for the function Pem, the method 
of maximising the a posteriori probability is equivalent to a semi-classical solution. The 
advantage of choosing a linear regulator is two fold. First, the search for the minimum 
is simply the solution of a system of linear equations; in non-linear minimisation it is no 
simple matter to correctly identify the global minimum Q]. Second, the linear problem 
is guaranteed to have a single minimum, whereas a general non-linear regulator may have 
multiple local minima, leading to complications analogous to the physics of phase transitions. 

Since previous work has demonstrated that for ^ T lattice computations match per- 
turbation theory P], we focus our attention on the region u < ttT. The linear relation 
between Gem and p^u means that we can assume p^M = P^em + ^Pemi where is the 
usual MEM prior in an ideal quark gas. At small uj this goes to zero faster than linearly in 
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co/tcT 



FIG. 2: /^pv{oj) obtained from fits to AGv (t) at T = 2Tc determined on a 12 x 26^ x 48 lattice 
with m/Tc = 0.03. Statistical errors obtained with = 32 are denoted by the bars, while the 
lines span thrice the range allowed by various systematic uncertainties as discussed in the text. 

uj and hence does not contribute to a/T [20]. By choosing to work with AGg^,^, this 
is removed from the problem, and we are freed to concentrate on the piece Ap^M, which 
contains all the information needed to extract cr. Then in eq. ((Tj) we use £ = 1, replace 
Gem by AG EM and Pem by Apem, and remove the term in p°^^. The upper limit of the 
integral was truncated to uj = 2mxT and the range divided into an uniform mesh of 
points. Varying n and A^^^ independently in the range 2 < n < 4 and 16 < A^^^ < 64 has no 
effect on the quality of the fit to the data (see Figure Q). 

Statistical errors on Apem are assigned by a bootstrap over the measured values of AGem- 
These are minor compared to uncertainties arising from algorithmic parameters. The latter 
are estimated by changing the integration method which is used to discretise eq. Q, the 
bin size in the integration, the integration limit, A'^, and the Bayesian prior specified by the 
matrix C. Thrice the combined uncertainty due to these four factors is shown as the band in 
Figure 121 within which ApEM lies. Even with this generous allowance for uncertainties there 
seems to be a peak in the spectral function at small u. For uj > tcT the spectral function 
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is roughly consistent with free field theory, but there is some evidence of a further peak at 
uj IGttT. The most important systematic uncertainty turns out to be related to control 
over the limit a 0. We found that the position of the peak and the slope at the origin 
change in going from A^^^ = 8 to 14. This phenomenon has been noticed earlier in the context 
of MEM |21j. A method which allows for better control of the continuum limit is required. 

For this we utilize a second class of Bayesian methods, in which the prior is a model of 
the observed bump in the soft part of the spectrum. Since Pem is real and non-singular for 
real uj, odd in uj, and non-negative for > 0, one can choose to work with the most general 
form which gives rise to a non- vanishing electrical conductivity, 

1 zY^^_ 7 



221. The constraint that 



where z = uj/T and with all 7„ and 6m real and non- negative 

^Pem — > at large uj is imposed by choosing M > N [2^ • We shall use the notation ( A^, M) 
to denote a particular choice of A^ and M. Bayesian techniques for parameter estimation 



24|. 



then proceed by choosing a priori probability distributions for each parameter 

The parameters in eq. (jH)) are most conveniently extracted by fitting to the Fourier trans- 
form of eq. (jni) over the Euclidean time t [2^ 

QemK^u, P,T)= f , (9) 

J llTX UJ — UJn 

where the Euclidean frequencies are uj^ = 2imTT, {1 < n < Nt on a lattice) and the 
path of integration over complex uj runs over the real line and is closed in the upper half- 
plane. The form of Pem in eq. (jH)) can then be used to express the Fourier coefficients in 
terms of the parameters, which can be determined either through a least squares method if 
1 + N + M < Nt/2, or a Bayesian method when there are more parameters than data. 

A particular simplification occurs for (0,M), since 70 is the only parameter that con- 
tributes for ujn = 0. In all these cases a/T = Cem{Xv ~ Xv')/37'^, where Xv = ^v(0,O,T) 



is the vector meson susceptibility [27[ obtained in QCD and the same quantity for 

an ideal quark gas on the same lattice [2^. Since the remaining parameters do not appear 
in this expression, their prior probabilities can be integrated out of the problem, without 
any assumptions about them. Such a marginalisation of the prior distribution is a general 
technique which has been demonstrated on other problems in the past 

The extraction using A^ = must be insufficient for T > lO^^T^, since it does not repro- 
duce the parametric dependence of a on g in weak-coupling theory, which is expected to 




o 


00 


m 


?^ 




00 




s 






Q. 


Q. 


00 

T3 




S 




II 

T3 




II 

O 





T/T 



FIG. 3: The electrical conductivity of the QCD plasma, cr, as a function of the temperature, T. 
The dimensionless quantity related to the soft photon emission rate shown on the right hand y-axis 
equals 6cr/T 20]. The bars denote statistical errors in the fit to the form in eq. Q with M = 2 
and N = 1. 



work at these temperatures. This can be improved by allowing for other values of {N,M). 
We have investigated the stability of our results by going to (1,2). Such a multiparameter 
fit moves the result up by 7%, which is within the statistical uncertainty. The electrical con- 
ductivity is thus reasonably stable, although it would be interesting in future to investigate 
its stability further, especially by using larger values of Nt. 

In principle, such an extraction of parameters other than 70 in eq. (jH)), allows us to proceed 
beyond the u = limit of the dilepton rate. As more parameters are determined, the shape 
of the soft dilepton spectrum is also better constrained. An interesting open question is of 
the number of Fourier coefficients needed to fix the shape of the dilepton spectrum. This 
question is related to the stability of the transport coefficient, and we plan a study in the 
near future to address this question. 

The estimates of a from the formula above are subject to lattice artifacts of order 
coming from Qy The continuum limit can then be obtained by an extrapolation in l/N^. 
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Finite volume effects turn out to be invisible within errors. Nor is there any visible quark 
mass dependence for small quark masses, since m/Tc = 0.03 and 0.1 give identical results 
within errors. We estimate a/T ^ 7Cem in the continuum limit of the temperature range 
we studied. Finally, have used the estimate of a to predict the soft photon emissivity of the 
QCD plasma in equilibrium, as shown in Figure 121 

In summary, we adopted a sequence of Bayesian techniques for the inverse problem of 
extracting spectral densities, Pem- In view of the results of we used the dispersion 
relations for Euclidean propagators after subtraction of the ideal gas values of the Euclidean 
temporal correlator, AGem- We observed that in the QCD plasma, in the temperature 
range 1.5 < T/T^ < 3, the spectral density is peaked at small energies. This peak was next 
analyzed using a parametrised form of the Bayesian prior and the electrical conductivity of 
the plasma was extracted. This was then used to predict the soft photon emissivity of the 
QCD plasma. 

Rough estimates of typical transport related time and length scales in the QCD plasma 
can be obtained by using the extracted value of the electrical conductivity in conjunction 
with the simple transport formula a = CEMnqTq/m {Cem is nothing but the average charge 
square: e^). If the number density of quarks, rig, is substituted by the corresponding entropy 
density, and the screening mass used for m, then one finds the quark mean free time ^ 0.3 
fm. This is also the time scale for the persistence of charge, isospin, strangeness and baryon 
number fluctuations in the plasma [5] . A possible experimental check could be to determine 
the mean free path of very soft off-shell photons {u ^ 0.15(27rT) ^ 200 MeV). These are 
only 20 times longer than Tq, i.e., about 6 fm. The fireball at RHIC may be marginally 
transparent to such photons, but at LHC the fireball size could be large enough to attenuate 
the intensity of very soft photons. Such an observation would constitute direct evidence for 
short mean free paths in the plasma. 

Some estimates of other transport coefficients can be obtained if one assumes that the 
mean free time of gluons is Tq/2, since they should be related by colour factors. Then simple 
transport formulae treated in the same approximation as before lead us to the estimate that 
the dimensionless ratio rj/S ^ 0.2, where S is the entropy density of the plasma and rj is 
the shear viscosity. This ratio is of the same order of magnitude as extracted from present 
heavy- ion data 29] and obeys a bound conjectured in j30|. It would be useful to make a 
direct measurement of the shear viscosity on the lattice. 
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Many interesting lines of research are relegated to the future. The dilepton emissivity 
away from = is a conceptually simple extension, but requires further numerical work, as 
explained before. Extending these measurements closer toward where correlation lengths 
grow larger ^] is of obvious importance, but outside the scope of this paper, as is the 
extension to dynamical QCD. The interesting question of the effectiveness of linear response 
theory, and hence of the Kubo formulae closer to Tc, can perhaps be probed using the 
non-linear susceptibilities defined in js^. 
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